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ABSTRACT 

The standard model of GRB afterglows assumes that relativistically 



O 
O 
O 

expanding material is decelerating due to interaction with the surrounding 



medium. The afterglows are well described by the synchrotron radiation from a 
forward shock, while the strong optical flash associated with GRB 990123 can 



O 

be attributed to the emission from a reverse shock. We give a detailed study 

on the reverse shock emission. The full light curves are calculated for a long 

and a short GRB cases. We discuss the lack of the prompt optical detections by 

2 ! ROTSE for GRB 981121 and GRB 981223. 
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O ' 1- Introduction 

The gamma-ray burst (GRB) afterglow is believed to involve a relativistically 

expanding fireball. The surrounding matter, which we will refer to as the ISM, influences 

the fireball shell after it has been collected enough and the considerable energy has been 

transferred from the shell to the ISM. The energy transfer is due to two shocks: a forward 

shock propagating into the ISM and a reverse shock propagating into the shell. 

The afterglow observations are fairly well described by the synchrotron emission from 
the ISM electrons accelerated by the forward shock, it is considered as a confirmation of 
the relativistic fireball. However, the current afterglow observations detect the radiation 
from several hours after the GRBs. At this stage, the Lorentz factor of the forward shock 
is not ultra-relativistic, less than ~ 10. Furthermore, the dynamics depends only on two 
parameters: the explosion energy and the ISM density. The afterglow observations provide 
neither verification of the extreme relativistic motion nor the properties of the fireball which 
can constrain models of the GRB source. 

The counterpart of the afterglow, the emission from a reverse shock was also predicted 
(Meszaros & Rees 1997; Sari & Piran 1999a). When a reverse shock crosses a shell, the 
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forward shocked ISM and the reverse shocked shell carry comparable amount of energy. 
However, the typical temperature of the shocked shell is lower since the mass density of 
the shell is higher. Consequently, the typical frequency from the shocked shell is lower. A 
prompt optical emission from GRB 990123 (Akerlof et al. 1999) can be regarded as this 
emission (Sari & Piran 1999b; Kobayashi & Sari 2000). 

The emission from the reverse shock is sensitive to the initial properties of the fireball. 
The observations can provide some important clues on the nature of the GRB source. 
Previous studies focused on the emission at the peak time. In this paper we calculate the 
full light curves for several frequency regimes. We in section 2 discuss the hydrodynamics of 
the relativistic fireballs on which the light curves highly depend. In section 3 we calculate 
the light curves for a long burst and a short burst case. We compare our estimates with 
the ROTSE observations in section 4. We estimate the initial parameters of the fireball of 
GRB 990123, and then we make some comments on the lack of the prompt detections by 
ROTSE for GRB 981121 and GRB 981223. In section 5 we give conclusions. 



2. Hydrodynamics of a Relativistic Shell 

Consider a relativistic shell with an energy E, a Lorentz factor rj and a width in 
laboratory frame Ao expanding into a surrounding medium (ISM) with a particle number 
density n\. When the shell sweeps a large volume of the ISM, it begins to be decelerated. 
The interaction between the shell and the ISM is described by two shocks: a forward shock 
propagating into the ISM and a reverse shock propagating into the shell. There are four 
region separated by the two shocks: the ISM(denoted by the subscript 1), the shocked 
ISM(2), the shocked shell material (3) and the unshocked shell material (4). Using the 
jump conditions for the shocks and the equality of pressure and velocity along the contact 
discontinuity, we can estimate the Lorentz factor 7, the pressure p and the number density 
n in the shocked regions as functions of three variables 711,714 and 77 (Blandford and McKee 
1976). 

There are two limits to get a simple analytic solution (Sari and Piran 1995). If the 
shell density is high 714 ^> 7/ 2 ni, the reverse shock is Newtonian which means that the 
Lorentz factor of the shocked shell material 73 is almost unity in frame of the unshocked 
shell material. It is too weak to slow down the shell effectively 73 ~ 77. On the other hand, 
if the density is low 774 <C i] 2 ni, the reverse shock is relativistic 73 ~ (711/774) ^(t^) 1 / 2 ^> 1 
and it considerably decelerates the shell material 73 ~ (n4/ni) 1 /*(r)/2) 1 / 2 <C 77. Once 
73 is determined, the density and the pressure in the shocked shell region are given by 
77 3 ~ (473 + 3)n 4 and p 3 ~ 47 2 77 1 m p c 2 /3. 
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In both cases, the time it takes for the reverse shock to cross a distance dx in the shell 
material can be given in a similar form up to a constant factor (Kobayashi, Piran & Sari 
1999), 

dR/c oc ^(ni/ni) 1 ' 2 dx/c, (1) 

where R is the radius of the shell. Since the motion of the shell is highly relativistic, we can 
regard R/c as time in laboratory frame. 

When a shell is ejected from a source, it has a high density n^ ^> i] 2 ni, so the reverse 
shock is initially Newtonian. However, as the shell expands, the shell density decreases. 
There is a possibility that the shock becomes relativistic during it is crossing the shell. Sari 
and Piran (1995) showed that, using the Sedov length I = (3E /AimirripC 2 ) 1 / 3 , if the shell is 
thick: A > l/2-rfl 3 , the reverse shock becomes relativistic at R N = l 3 ^ 2 /A ' rj 2 before the 
crossing at Ra = / 3//4 A . If the shell is thin: A < l/2rj 8 ^ 3 , it is dense enough to keep the 
reverse shock Newtonian. The shock becomes mildly relativistic only when it just crosses 
the shell at i? 7 = l/rf' 3 . 

The typical burst energy is about 10 52 ergs. The ISM density has a typical values of 
5 protons/cm 3 . The thick solid line in figure [I] separates the thick shell case (upper right) 
and the thin shell case (lower left). According to the internal shocks model the duration of 
a GRB T is given by the shell width A /c, the thick shell cases correspond to relatively 
long bursts. 



2.1. the Thick Shell Case 

In a thick shell case the reverse shock becomes relativistic at Rn and it begins 
to decelerate the shell material. Since the shell density decreases with radius as 
n 4 ~ ra 1 / 3 /r/ 2 A -R 2 , using equation [[], one finds that the number of shocked electrons N e is 
proportional to R 2 . The scalings of the hydrodynamic variables in terms of the observer 
time t = R/2cy 2 are 

73 ~ (//A ) 3/8 (4ct/A r 1/4 ; n 3 ~ 8 1 3 n 1 /r ] oc t~ 3 '\ (2) 

p 3 ~ 4 73 2 nim p c 2 /3 oc r 1/2 ; N e ~ N ct/A (3) 

where N = E/r]m p c 2 is the total number of electrons in the shell. Note that the Lorentz 
factor 73 does not depend on its initial value r] after the reverse shock becomes relativistic. 

After the shock crosses the shell at t — Aq/c, the profile of the forward shocked ISM 
begins to approach the Blandford-McKee (BM) solution (Kobayashi, Piran & Sari 1999). 
Since the shocked shell is located not too far behind the forward shock, it roughly fits the 
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BM solution. The author and Sari (1999) numerically showed that the evolution is well 
approximated by the BM solution if the relativistic reverse shock can heat the shell to a 
relativistic temperature. Then, we apply the BM scalings to the evolution of the shocked 
shell. The number of the shocked electrons is constant after the shock crossing because no 
electron is newly shocked. 

73 oc r 7/m , n 3 oc r 13/16 , p 3 oc r n/12 , N = constant. (4) 



2.2. the Thin Shell Case 

In a thin shell case the reverse shock is too weak to decelerate the shell effectively, the 
Lorentz factor of the shocked shell material is almost constant during the shock propagates 
in the shell. Due to a slight difference of the velocity inside the shell, the shell begins to 
spread as A ~ R/rf around R = i] 2 Ao. Then, the density ratio decreases as n^/rii ~ (l/R) 3 . 
The scalings before the shock crosses the whole shell at t 7 = l/2cr] 8 ^ 3 are given by 

73 ~ V, n 3 ~ 7nm 2 (t/t 7 )- 3 : p 3 ~ 4^n imp c 2 /3, N e ~ iV (£/t 7 ) 3/2 . (5) 

In the above thick shell case the spreading effect was not important since it happens after 
the shock crossing. 

The Newtonian reverse shock can not heat the shell to a relativistic temperature which 
the BM solution assumes, then, we are not able to use the BM solution. However, we 
derived scaling laws for a cold shocked shell, assuming that the Lorentz factor is described 

4/3 

by a power law ex R~ 9 and that the shell expands adiabatically p 3 oc n 3 with a sound 



speed ~ slPs/n^ in the shell's comoving frame. We numerically showed that the scalings 
with g ~ 2 fit the evolution (Kobayashi & Sari 2000). 



73 oc 



r 2/5 , n 3 oc r 6/7 , p 3 oc t~ 8/7 , N = constant. (6) 



3. Light Curves of the Reverse Shock Emission 

We consider now the synchrotron emission from a reverse shocked shell. The 
shock accelerates electrons in the shell material into a power law distribution: 
N^ejdje oc %^d r y e (7e > 7 m ). Assuming that a constant fraction e e and e B of the 
internal energy go into the electrons and the magnetic field respectively, one finds that the 
typical random Lorentz factor of the electrons and the magnetic field evolve as 7 m oc e e p 3 /n 3 
and B 2 oc esp 3 . 
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The spectrum is given by the broken power laws discussed in Sari, Piran & Narayan 
(1998). In this paper we neglect the self absorption since it does not affect the optical 
radiation in which we are interested. Then, it has two breaks at the typical synchrotron 
frequency v m oc .6737^ and at the cooling frequency v c oc 1/ 'B 3 ^^ 2 which is the synchrotron 
frequency of electrons that cool on the dynamical time of the system. The peak flux is 
obtained at the lower of the two frequencies. Let N e and D be the total number of the 
shocked electrons and the distance to the observer respectively. The observed peak flux 
density evolves as F u ^ max oc N e Bj 3 /D 2 . 

The spectra do not depend on the hydrodynamics of the shell. However, the light curve 
at a fixed frequency depends on the temporal evolution of the break frequencies v m and v c 
and the peak power F u ^ max . These depend on how 73, n 3 ,p 3 and N e scale as a function of 
t. We will apply the adiabatic evolution discussed in section 2 to the shell evolution. It is 
justified if the fraction of the energy going to the electron is small e e C 1 or if we are in the 
regime of slow cooling u m < v c where the electrons forming the bulk of the population do 
not cool. 



3.1. the Thick Shell Case 

A reverse shock crosses a thick shell at ~ A /c, the peak time of the emission from the 
reverse shock is comparable to the GRB duration T. Using the estimates in section 2 we 
obtain the break frequencies and the peak flux at the shock crossing time, 

u m ~ 7.3 x 10 14 e^ 2 n^ 3 2 00 Hz, (7) 

v c ~ 9.4 x 10 15 e~ B l /2 E^ /2 n^Tj 2 Hz, (8) 

F no r>-2 V 2 Z7> 5 / 4 V 4 -1t-i-3/4 T / n x 

fu,max ~ U.d JJ 28 e m h b2 n 15 77 300 i 100 Jy, (9 J 

where we have scaled the parameters as E 52 = -E/10 52 ergs, e e0 = e e /0.6, e^o = es/0.01, 
^1,5 = n.i/5 cm -3 , 77300 = r//300, T100 = T/100 sec and D 2 % = D/10 28 cm. The equipartition 
values e e = 0.6, e^ = 0.01 and n\ = 5 protons/cm 3 are inferred for GRB 970508 (Wijers & 
Galama, 1999; Granot, Piran & Sari 1999). The scalings before and after the shock crossing 
are given by 

t <T : v m — constant, u c oc t~ l , F umax oc t 1 ' 2 , (10) 

t > T : z/ m oc t" 73 / 48 , v c oc t 1 / 16 , F Utmax oc t" 47 / 48 . (11) 

It is interesting that v m is constant during the shock crossing. 

The spectrum is the slow cooling throughout the evolution if v m < u c at T. The shaded 
region in figure [1| shows the corresponding parameter region. The flux at a given frequency 



v evolves as 
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f t 1/2 v < v 

FJt<T) oc c (12) 

I constant v > v c 



'in 



£- 17 / 36 V <v n 

I'M- V) x <Jt-(^+ 21 )/ 96 v m <v<v cut (13) 

v > u cut 

The flux at a frequency above v c disappears at T because no electron is shocked anymore. 
This cut off frequency v cut decreases as z/ c (T)(£/T)~ 73 / 48 due to the adiabatic expansion of 
the fluid. The index — (73p + 21)/96 is ~ —2.1 for the standard choice p = 2.5. 

The light curves are different among three frequency regimes separated by the two 
frequency, v m and v c at T. The typical light curves are shown in figure ||a. The flux 
initially increases at all frequencies as t 1 / 2 , but the high frequency light curve flats when the 
cooling frequency crosses the given frequency. At T, the flux begins to decay as £~ 17 / 36 or 
£-(73?5+2i)/96^ pj na ]iy it vanishes when v cut crosses the given frequency v at 

, 7nn -72/73 ,7-24/73 -48/73^49/73/ v \ -48/73 

t ~ 700 t m E 52 n lj5 Tioo 7 ^ - - — 4 — j sec. (14) 

With a parameter set in the upper right region of the dashed line in figure [3], the 
spectrum changes from the slow cooling to the fast cooling during the shock crossing. The 
light curves are different among three frequency regimes separated by v m and u c at T. The 
typical light curves are shown in figure |2|b. We neglected the initial slow cooling phase in 
the figure since the transition happens at very early time. The flux at a given frequency 
above v c is constant and below u c it increases as £ 5 / 6 . When the shock crosses the shell, 
above v c it vanishes. After that, below v c it decreases as £~ 17 / 36 until v cut crosses the given 
frequency. 



3.2. the Thin Shell Case 

As a thin shell case, the shell width Aq/c should be smaller than t 7 . Then, separation 
is expected between the GRB and the peak of the emission. The corresponding parameter 
region is lower left of the solid line in figure [l|. The break frequencies and the peak flux at 

, o ,-,1/3 -1/3 -8/3 . , 

t 7 ~ 3 E^ n 1)5 ' 77300' sec are given by 

u m {ty) - 9.6 x 10 14 e^V^oo Hz, (15) 

i/efo) - 4.0 x 10 16 e^V 3 n^g Hz, (16) 

F U) max{t~,) ~ 5.2 D^e 1 ^E 52 n l 1 / 5'n m o Jy. (17) 
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Note that these do not depend on the initial shell width A itself though A /c should be 
smaller than £ 7 . The scalings before and after £ 7 are as follows, 

t < t 7 : z/ m oc t 6 , v c oc £~ 2 , F„ imoa; oc £ 3/2 , (18) 

t > t 7 : z/ m oc t" 54 / 35 , ^ c oc t 4 / 35 , F„, maa; oc t" 34 / 35 . (19) 

If the spectrum is the slow cooling throughout the evolution, the behavior of the light 
curves are different among the three frequency regimes separated by v m and v c at t 7 . Figure 
|3]a shows typical light curves. The flux initially increases in all regimes very rapidly as 
£3p-3/2^ ^6 f or p _ 2.5). The slope changes when v m or v c crosses the given frequency. 
After t 7 , the flux decays as t^ 16 / 35 or t~( 27 P+ 7 )/ 35 (~ t~ 2A for p = 2.5). The emission at a 
frequency z/ < z/ c (t 7 ) disappears when the cut off frequency crosses it at 

* ~ 50 e-£l »E^*Zr / *%X m ( 5 - ^ H J " 35/M »c (20) 

If z/ m is higher than v c at t 7 , the spectrum changes from the slow cooling to the 
fast cooling during the shock crossing. However, it requires a large Lorentz factor 
r) > 8 x 10 4 ^eo e Bo^52 n i^- Then, the transition happens only if the shell is extremely 
thin (see figure [TJ), then it is hard to detect such a prompt emission. We show the 
possible light curves for completeness. The behavior is different among the following four 
frequency regimes: v < u c (t 7 ), ^ c (^ 7 ) < v < Uo, v < v < u m (t 7 ) and v > u m (t y ) where 
^o = ^m(^o) — v c (to) and to is the transition time from the slow cooling to the fast cooling. 
The typical light curves are shown in figure |3|b in which the very early slow cooling phase is 
neglected. 



4. Observations 

We in this section compare our estimates with the ROTSE observations. First, we 
determine the parameters of the fireball of GRB 990123 from the observations. Then, we 
make some comments on the lack of the prompt optical emission from GRB 981121 and 
GRB 991223 in the context of our model. 



4.1. GRB 990123 

GRB 990123 is a very bright burst which fluence is about 100 times that of a median 
BATSE burst. Absorption lines in the optical afterglow gives a lower limit of the redshift 



z > 1.6, the energy required to produce the bright GRB is enormous, ~ 10 54 ergs for an 
isotropic emission. ROTSE detected a strong optical flash during the burst (Akerlof et al. 
1999). About 50 sec after the GRB trigger, it reached to a peak of ~ 1 Jy and then decayed 
with a slope of a power law index ~ 2. 

The fireball of GRB 990123 is likely to be a marginal case (Kobayashi & Sari 2000). 
However, such a marginal case behaves very much like the thin shell case, since the shell 
heated by a mildly relativistic shock becomes cold soon. The observed decay is in good 
agreement with the theory t~ 2A . This means that the relation among u m , v c and the 
observed optical frequency ^~5x 10 14 Hz is v m < vr < v c at the peak time. Furthermore, 
the strong optical emission implies that vr is close to v m (Sari & Piran 1999b). 

It is well-known that the peak time t pea k is sensitive to the Lorentz factor of the 
shell, so we obtain 73 ~ 180 E 5i n l J (t pea k/50 sec)~ 3//8 . Since the reverse shock is mildly 
relativistic, the Lorentz factor at the peak time should be close to the initial Lorentz factor 



which can be estimated by using equation 20. The last detection of the optical flash by 



ROTSE was ~ 600 sec after the GRB trigger. Using the equality of the Lorentz factors, we 
obtains 77 ~ 270 and m ~ 0.2 protons/cm 3 . With these parameters the typical synchrotron 
frequency at the peak time is ~ 1.6 x 10 14 Hz which is close to vr and it is consistent with 
the initial assumption. 

Though we can not give a strong argument on the raising part of the light curve since 
the observation is sparse, the power law index of the slope calculated from the first two 
ROTSE data is 3.4. The sparseness can make only the index smaller, then the real index is 
at least larger than that of a thick shell. 



4.2. GRB 981121 and GRB 981223 

The theory succeeded for GRB 990123, but the optical flash was detected only for 
it so far. Akerlof et al. (2000) reported no detections of the optical flashes to six GRBs 
with localization errors of 1 deg 2 or better. Especially, GRB 981121 and GRB 981223 
are the most sensitive bursts in the sample. If the optical flashes are correlated with the 
GRB fluences, the optical emission should be more than 2 mag over the ROTSE detection 
thresholds. The thick lines in figure |5|a and [|a show the light curves expected from the 
GRB fluences, /auBQsim = 7.0 x 10 -6 ergs/cm 2 , f G RB98i223 = 1.3 x 10" 5 ergs/cm 2 and 
fcRB990i23 = 1.0 x 10~ 4 ergs/cm 2 (Akerlof et al. 2000). 

GRBs are produced by the internal shocks, while the afterglows and the optical flashes 
are due to the external shocks. The lack of a direct scaling between the GRB and the 
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afterglow is an evidence of the internal shocks model. However, the energy radiated in the 
GRB phase E 7 is still correlated with the blast wave energy E estimated from the x-ray 
afterglow (Freedman and Waxman 1999). According to their analysis, the ratios E 7 /E 
of 13 events are in between 0.1 and 6, while E 1 ranges from 3 x 10 51 ergs to 10 54 ergs. 
Furthermore, E^/E ~ 3 of GRB 990123 is a relatively large value in the sample. Therefore, 
the scaling with the GRB fluence does not overestimate D~ 2 E in general. 

GRB 990123 was an exceptionally energetic burst. The energies of GRB 981121 and 
GRB 981223 should be considerably lower to explain the lower fluences, otherwise the 
sources are extremely far, z=10 for GRB 981121 and z=6.5 for GRB 981223 assuming 
h = 0.65, f^o = 1 and Ao = 0. Hereafter we assume E = 10 52 ergs for the two bursts. Since 
GRB 990123 is a marginal case, a burst with a lower E and a comparable ni, rj and T 
is classified into the thick shell case in which F v>max is proportional to E b / A instead of E. 
Then, the peak flux for E = 10 52 ergs is lower by a factor of ~ 3 than that just scaled by 
the GRB fluence. The thin lines in figure |5|a and ||a depict the corrected light curves which 
are still above the thresholds. The durations of GRB 981121 and GRB 981223 are 54 sec 
and 30 sec respectively and comparable to that of GRB 990123. 

A possible solution to the problem is to assume that the reverse shock energy is 
radiated at a non-optical frequency, v m <C Vr or v m ^$> vr. The typical frequency v m is 
proportional to e 2 e B n/ rj 2 , but the values of e e and e# are determined by the micro-physics 
and are likely to be universal. The difference of v m should be due to that n-y and rj of the 
bursts are different from the "canonical" values, n\ = 0.2 protons/cm 3 and rj = 270. 

If n\ and rj are smaller than the canonical values, v m is lower than vr since v m ~ Vr 
for GRB 990123. Using the normalization by the GRB fluence, the peak flux is 



i<VR< 



, c (T)~1.0/o^ /4 T - 3/4 < /4 ^ 2 Jy (21) 



where the subscript denotes that the parameters are scaled by the values of GRB 990123, 
/o = //10~ 4 ergs cm -2 , E = .E/10 54 ergs, T = T/50 sec, n = ni/0.2 protons cm -3 and 
Vo — 77/270. The fireball is classified into the thin shell case if we assume very small ri\ or 



(p+l)/4 p 



77, thus the dependence of the peak flux on n\ and rj changes to F VR {t^) oc Uq rj^. F t 



I'R 



and F vmax at the peak time are plotted in figure [| as functions of rj and n\ in the case 
of GRB 981121. The no detections by ROTSE give upper limits on n\ and rj. Assuming 
m = 0.2 protons/cm 3 , we get rj < 135 for GRB 981121 and rj < 65 for GRB 981223, or 
assuming rj = 270, n\ < 0.07 protons/cm 3 for GRB 981121 and rix < 0.02 protons/cm 3 for 
GRB 981223. The light curves with a low n\ or rj are shown as the thin lines in figure ||b,c 
and ||b,c. 
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If we assume a larger n\ or 77 than the canonical values, the peak flux is given by 

Fu R <u m (T) ~ 1.0 f Q El'V A nT^ Jy. (22) 

The spectrum changes to the fast cooling if n\ or rj is very large, the dependence on m 
and 7] changes to F VR oc n ' rj . For a further larger %, z/ c can be lower than z/r, then 
the dependence becomes n . The no detections give lower limits, assuming n\ — 0.2 
protons/cm 3 , r? > 400 for GRB 981121 and r] > 600 for GRB 981223, or assuming r] = 270, 
ni = 2x 10 5 protons/cm 3 for GRB 981121 and m = 4 x 10 6 protons/cm 3 for GRB 981223. 
The light curves for a large n\ or r\ are shown as the thick lines in figure [|b,c and |]b,c. 

A large ISM density is an unlikely reason to explain the no detections, because we 
need to require several order larger density than the canonical and with the large density 
the peak power F v>max itself is very large (see figure f|b). However, there is a possibility of 
the extinction. Though we have normalized the optical flux according to the GRB fluence, 
gamma-rays do not suffer any kind of extinction. A half of x-ray afterglow bursts do not 
have optical afterglows, it might be due to the extinction. The reverse shock radiates at a 
closer region to the inner engine, it might be crucial. Within the six bursts to which the 
ROTSE group reported the no detections, only the location of GRB 980329 was determined 
precisely by BeppoSAX, and the optical afterglows were observed hours later. However, the 
ROTSE observation on this burst was not so sensitive. Future observations will provide 
some information to this effect. 



5. Conclusions 

We have constructed the full light curves of the reverse shock emission for a short burst 
(thin shell) and a long burst (thick shell). The typical synchrotron frequency increases 
rapidly as t 6 in the thin shell case while it is constant in the thick shell case. For a plausible 
moderate Lorentz factor, upto a few thousand, the synchrotron spectrum is the slow cooling 
throughout the evolution, with which spectrum we find that the flux must rise initially 
steeply as £ 3 p -5 / 2 or £ 3 p~ 3 / 2 in the thin shell case, and slowly as t 1//2 in the thick shell case. 

The rise ends when the reverse shock crosses the shell. Only the exception is the low 
frequency emission from a thin shell which already begins to decrease when the rapidly 
changing typical frequency crosses the observe one. In the decay phase, the light curves are 
similar for both cases, though the hydrodynamics are very different. Then, the detection 
before the peak should be more useful to give a constraint on the initial properties of the 
fireball. 

The prompt optical emission from GRB 990123 is well described by the reverse shock 
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emission. The observations enabled us to determine the initial Lorentz factor and the ISM 
density. We found that GRB 990123 is a "luck" burst. Besides the exceptionally large 
energy, it has an optimized Lorentz factor to produce a bright optical flash. As the initial 
Lorentz factor increases, the peak power, F v ^ max oc 73/77, initially rises. However, with a 
moderate initial Lorentz factor the evolution changes to the thick shell case in which the 
Lorentz factor at the peak time no longer depends on the initial Lorentz factor. Then, 
the peak power drops since the number of the electrons in the shell continues to decrease. 
Therefore, the marginal case gives the brightest emission. GRB 990123 is the marginal case 
and the typical frequency v m oc n/ rf just comes to the optical band. 

The lack of the prompt optical detections by ROTSE for GRB 981121 and GRB 981223 
does not give strong constraints on the initial Lorentz factors or the ISM densities. If the 
Lorentz factor is slightly different from that of GRB 990123, the peak flux becomes lower 
than the ROTSE thresholds. It is also possible to explain the no detections by a lower ISM 
density. 

The author thanks Re'em Sari for helpful discussions, and Robert Kehoe for providing 
with the ROTSE data. This work was supported by the Japan Society for the Promotion 
of Science. 
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Fig. 1. — Allowed parameter region: The thick solid line separates the thick shell case (upper 
right) and the thin shell case (lower left). The slow cooling region for the thick shell case 
is lower left of the dashed line. That for the thin shell is 7] < 8 x 10 4 . The peak time of 
the reverse shock emission is t pea k = max[t 7 ,T], the thin solid line depicts t pea k = 1 sec. 
E = 10 52 ergs, n\ — 5 protons/cm 3 , e e = 0.6 and et = 0.01 are assumed. 



-14- 



1(T 




10 



"3 



10 



10 L 



2 










b 






v= 


c 


I 


const. 




t - 17/36 




y^ t 5/6 










- 


4 








const. 




- 












:T 



10' 



10' 

t(sec) 



1(T 



10 H 



Fig. 2.— Light curves: thick shell case. (a) Slow cooling case: i] = 300. thin 
solid (i/ = 10 13 Hz< i/ m (T)), thick solid (i/ m (T) < v = 10 15 Hz< u c (T)) and dashed 
[y = 10 17 Hz> u c (T)). (b) Fast cooling case : r] = 10 4 . thin solid [y = 10 15 Hz< u c (T)), 
thick solid(z/ c (T) < v — 10 17 Hz< u m (T)) and dashed [y = 10 19 Hz> u m (T)). v = v m and 
v = v c show the times when the observe frequency v is crossed by v m and v c respectively. 
E = 10 52 ergs, n\ — 5 protons/cm 3 , e e = 0.6, ej, = 0.01 and T = 100 sec are assumed. 
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Fig. 3.— Light curve: thin shell case. (a) Slow cooling case: rj = 300. Thin 



solid(z/ = 10 Hz < u m (tj)), thick solid(z/ m (t. 



< v 



10 5 Hz < jy c (t-y)) and dashed 



[y = 10 17 Hz > z/ c (£ 7 )). (b) Fast cooling case, thin solid (v < udt^)) , thin solid 
(i/ c (t 7 ) < v < z/ ), thick dashed (u < u < ^ m (t 7 )) and thin dashed {y > z/ m (t 7 )). v = u m 
and v = v c show the time when the observed frequency is crossed by v m and v c respectively. 
E = 10 52 ergs, ri\ = 5 protons/cm 3 , e e = 0.6 and tb = 0.01 are assumed. 
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Fig. 4. — GRB 981121: the optical flux F VR (thick) and the peak power F v ^ max (thin) at the 
peak time max[£ 7 ,T]. v m = vr shows the points for the canonical parameters. 
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Fig. 5.— GRB 981121: the ROTSE detection thresholds (segments) and the theoretical 
light curves, (a) E = 10 52 ergs (thin) and E = 10 54 ergs (thick), (b) i] = 100 (thin) and 
f] = 1000 (thick), (c) rii = 10~ 3 protons/cm 3 (thin) and n\ = 10 6 protons/cm 3 (thick). 
rii = 0.2 protons/cm 3 , r] = 270, E = 10 52 ergs and T = 54 sec are used if the values are not 
specified. The thresholds are calculated assuming that the ROTSE magnitude corresponds 
to that in the V band. 



-18- 




E 



1(T 




10' 



10' 

t(sec) 



1(T 



Fig. 6. — GRB 981223: the ROTSE detection thresholds (segments) and the theoretical light 
curves, (a) E = 10 52 ergs (thin) and E = 10 54 ergs (thick), (b) rj = 60 (thin) and 77 = 1000 
(thick), (c) rii = 10 -3 protons/cm 3 (thin) and n\ = 10 8 protons/cm 3 (thick). n\ = 0.2 
protons/cm 3 , rj = 270, E = 10 52 ergs and T = 30 sec are used if the values are not specified. 



